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1 ) Label root voxel with integer 0* 

2) Construct a queue and line up the root in the queue, 

3) If (There is at least one voxel in the queue) 

Serve the voxel x on the top of the queue: 
For (each of x's 26-connected neighbor voxel y) { 
If (y in the volume and has not been labeled yet) { 
Line up the y in the queue* 

/* label the voxel y */ 
Set dist =999999} 

For (each of y's 26-connected neighbor voxel z) { 

If (z in the volume and has been already labeled with an 
integer of n z ) { 

d z =n z +d (y,z)i 

where d (y,z) is 10, 14 or 17 if the Euclidean distance 
between y and z is 

1,^2, or ^3, 

respectively* 

If {dist>d z ) { 

label y with integer dist^ 
dist = d z * 

} 

} 

} 

} 

^ x leaves the quere; 
Else { 

end of calculating the distance map. 
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